Comparing non-perturbative models of the breakup of neutron-halo nuclei 
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Breakup reactions of loosely-bound nuclei are often used to extract structure and/or astrophysical 
information. Here we compare three non-perturbative reaction theories often used when analyzing 
breakup experiments, namely the continuum discretized coupled channel model, the time-dependent 
approach relying on a semiclassical approximation, and the dynamical eikonal approximation. Our 
test case consists of the breakup of ^^C on Pb at 68 MeV/nucleon and 20 MeV/nucleon. 
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I. INTRODUCTION 



Due to the proximity to the particle threshold, loosely- 
bound nuclei dissociate easily during collisions with 
nuclear targets. Consequently, they are often stud- 
ied through breakup reactions, in which the loosely- 
bound particle (s) dissociates from the core of the nucleus 
through interaction with the target. In the following, we 
focus on elastic breakup, i.e., a reaction in which the tar- 
get is left in its ground state and all projectile fragments 
are detected in coincidence after dissociation. The use of 
breakup reactions for extracting properties of exotic nu- 
clei is numerous and varied, including one-neutron halo 
systems [lH3|, configuration- mixed systems [J, |5|, two- 
neutron halo systems [a-Qj as well as proton rich sys- 
tems pj|4ll| . While the shape of the energy distribution 
can tell us about the separation energy and the angu- 
lar momentum of the ground state, the magnitude of the 
cross section is related to the asymptotic normalization 
of the ground state [l^l- In addition, for two-particle 
halo systems, one also expects to obtain information on 
the correlations in the valence pair [13| . More recently, 
breakup reactions have proven to be a useful tool in ex- 
ploring nuclei beyond the dripline and studying decay 
modes of resonant states J14| - [l6j . 

If the reaction is dominated by the electromagnetic in- 
teraction, it is possible to connect the breakup cross sec- 
tion with the capture cross section [13, [3 • This method, 
known as the Coulomb-dissociation method, is of interest 
to astrophysics because it can provide radiative-capture 
cross sections at very low relative energies where a di- 
rect measurement is not feasible. It also gives access to 
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neutron-capture cross sections by unstable species, which 
are impossible to measure in the laboratory. It has been 
applied to a number of cases 0, m, [l8l - l20 |. Providing 
confidence that the Coulomb dissociation method works, 
neutron capture cross sections for ^'^C{n,j)^^G were ex- 
tracted from the Coulomb dissociation data [3| using two 
independent methods [2l|, I24I and perfect agreement was 
obtained when compared to direct measurements [23|. 
Similar efforts have been performed for the breakup of 

The common feature of all the above mentioned ex- 
periments is their need for a reliable reaction model in 
the analysis. One needs to be careful with separating 
nuclear and Coulomb processes, often nuclear- Coulomb 
interference is important and the dynamical effects in the 
continuum are crucial [18, 2i, 28 - 31] . These results im- 
ply that, in general, perturbative approaches are not ac- 
curate enough for a reliable analysis of Coulomb-breakup 
measurements. 

Hand-in-hand with the experimental advances, a num- 
ber of non-perturbative breakup theories have been de- 
veloped improving the method by which breakup re- 
actions are studied (see Ref. [32| for a recent review). 
The many new developments rely on different approx- 
imations, have separate advantages and shortcomings, 
and vary also in the level of complexity. We believe it 
is timely to compare these theories and understand the 
level of accuracy of the approximations made. In this 
work we compare the most common non-perturbative ap- 
proaches to describe the breakup of a one-neutron halo 
nucleus that can be approximated by a two-body cluster. 
These are: i) the continuum discretized coupled channel 
method (CDCC) [3J,l3J], which is fully quantal and does 
not make approximations in the projectile-target dynam- 
ics, ii) the time-dependent approach (TD) l35[-| 38i . which 
is based on a semiclassical approximation [39| that de- 
scribes the projectile-target relative motion by a classical 
trajectory, and iii) the dynamical eikonal approximation 



(PEA) [401 141| which rcHes on the cikonal approximation 

All three theories are built on the same three-body de- 
scription of the reaction: the projectile P, described as a 
valence neutron / loosely bound to a core c, impinges on 
a target T considered as inert. The effective interaction 
between c and / is adjusted to reproduce known proper- 
ties of the projectile, while the interactions between the 
projectile fragments and the target are simulated by opti- 
cal potentials fitted to elastic-scattering data for the c-T 
and f-T systems. 

For the type of reactions we are interested in here, 
the CDCC method is the most accurate method avail- 
able on the market but it is also the most computation- 
ally intensive and requires elaborate model-space stud- 
ies. Developments that go beyond the inert-core and/or 
inert-target approaches, or extensions to A^-body projec- 
tile clusters {N > 2) are compromised by computational 
limitations. On the opposite, TD and DEA are not com- 
putationally intensive and are rather straightforward to 
set up. The question is whether these approximations 
can do a good job for the reactions of interest. To answer 
this question one needs to quantify the level of accuracy 
of the approximations introduced. 

For a meaningful comparison, it is necessary that all 
three methods start from the exact same three-body 
Hamiltonian. Typical breakup observables are then com- 
pared to quantify the accuracy of the various approxima- 
tions. The test case chosen is 208pb(i5c, "C n)'^°^Ph, 
a case where we expect the three-body description to 
be adequate. Our study is performed at two energy 
regimes, one at a typical energy in fragmentation facil- 
ities (68 MeV/nucleon) Q and the other at the higher 
energy limit of ISOL facilities (20 MeV/nucleon). We 
ignore for practical reasons the effect of relativity. 

In Sec. In] we briefly summarize the three methods un- 
der scrutiny. The model inputs are given in Sec. IIIIl In 
Sec. lIVI the results for breakup are presented, and conclu- 
sions are drawn in Sec. |Vl The details about the calcula- 
tions in all three models can be found in the Addendum 
provided as supplemental material. 



II. BRIEF THEORETICAL DESCRIPTION 

A. Common frame w^ork 

To study the breakup of a projectile P into a core c 
and a valence neutron / on a target T, we start from the 
(non-relativistic) three-body Hamiltonian 

H3b{R, r)=fjf^ + Ho{r) + Uct{Rc) + UfriRf), (1) 

expressed in the set of coordinates illustrated in Fig. [T] 
In Eq. ([T|), Tjj is the kinetic-energy operator for the P-T 
relative motion. The two-body Hamiltonian Hq describes 
the internal structure of the projectile 



where TV is the c-f kinetic-energy operator and Vcf is an 
effective potential, modeling the c-f interaction. This po- 
tential is adjusted to reproduce the bound-state spectrum 
and low energy scattering states of the projectile. The 
optical potential Uct (U/t) describes the elastic scat- 
tering of the core (valence neutron) by the target and 
contains a Coulomb part and a nuclear part. 




FIG. 1: Set of coordinates used in the reaction modeling. The 
longitudinal Z and transverse b components of R are shown 
as well. 

In all three methods, a partial wave expansion for the 
projectile states is used: 



IjIM 



4>r\r) 



(r) 



[Yiir)®X,]^(SXi^ 



IM 



(3) 



where K is a spherical harmonic [43| and Xs are spinors. 
The quantum number I is the orbital angular momen- 
tum of / relative to c, s{Ic) is the spin of the fragment 
/(c), and the total angular momentum of the projectile 
is / with projection M. We denote by 4>o the projectile 
bound state of (negative) energy Eq. For simplicity in 
this formulation we consider only one bound state. In 
this manner, all the other eigenstates of Hq correspond 
to positive energies E = h /2fXcfk^, with ficf the c-f 
reduced mass. They describe the c-f continuum. Of 
course, the formulation can be easily extended to include 
bound excited states and we do include these in the ap- 
plication presented in Section [TVl 

Within this framework, the study of the P-T collision 
reduces to solving the Schrodinger equation 



i/3b*(fl,r) = ^tot*(i2,r) 
with initial boundary condition: 

^{R,r) ^ e'^'°^0o(r). 



(4) 



(5) 



Hoir)^fr + Vcf{r), 



(2) 



where the initial P-T momentum TiKq is assumed along 
the Z axis. Its norm is related to the total energy 
Etot = ^^KQ/2fipT + Eq, with fipT the P-T reduced 



mass. There are different assumptions used in the treat- 
ment of the full three-body wave function *!'. We capture 
the essential features in the following subsections. 



B. Continuum discretized coupled channel method 

The full three body wave function can always be ex- 
panded in terms of the complete set of projectile states 
as: 
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such that the momentum hk between the internal mo- 
tion of c -f- / is related to the momentum ?iK between 
the projectile center of mass and the target through en- 
ergy conservation E'tot = h^K'^/2^PT + h^k'^/2fj,cf- An 
expansion involving an integral over momentum is not 
tractable, so in CDCC a discretization of the projectile 
continuum is performed {33i . i44] . There are various ways 
of performing this discretization, and here we will use the 
so-called average method whereby (f)^ is replaced by 

its average over a momentum bin [fcp_i,A;p], ^i-'^*^ |44| . 
In this method the three-body wave function is approxi- 
mated by 
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with p ~ corresponding to the initial ground state and 
p > I corresponding to the bin wave functions. The 
sum runs up to N, which is associated to the maximum 
projectile excitation energy i?max considered in the model 
space. In the end, the method needs to be independent 
of discretization and model space, and thus ii^max needs 
to be large enough and the bin width needs to be small 
enough to accurately describe the process of interest [4J] . 
When introducing expansion ([7]) into the full three- 
body equation ([¥]), and after integrating over the angu- 
lar variables and r. one arrives at the following coupled- 
channel equations in R |4J] : 
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where L is the P-T relative angular momentum, Ep is 
the midpoint energy of bin p, and a is the index for the 



channel {pljIL}. 
defined by, 

v±nR) - 



The coupling potentials Vj^^°' (R) are 



where Jtot is the total angular momentum resulting from 
the coupling of / and L. Equation ([5]) is solved with 
scattering boundary conditions at large distances: 
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where a^ is the entrance channel, and H are Coulomb 
Hankel functions [ijl ■ Breakup observables are then cal- 
culated from the resulting S matrix [SJ, 14^. In the 
present work, we use the code FRESCO to numerically 
solve the set of coupled equations (|8]) [43 ■ The parame- 
ters of our calculations are given in the Addendum pro- 
vided as a supplemental material of this article. 



C. Time-dependent model 

It can be very demanding to solve the coupled-channel 
equations ([5|) numerically. To reduce the computational 
cost, other models have been developed. In the semiclas- 
sical approximation, the P-T relative motion is approxi- 
mated by a classical trajectory R{t) |35l - (38t . Along that 
trajectory, the projectile experiences a time-dependent 
potential that simulates its interaction with the target. 
Assuming a quantal description of the internal structure 
of the projectile, this approximation leads to the time- 
dependent Schrodinger equation. 



d 



^^f^^'^'it, b, r) = [Ho + VpT{t, r)] vl/^"(i, 6, r), (11) 

where b is the impact parameter characterizing the tra- 
jectory. The time-dependent potential VpT appearing in 
this equation is the sum of the optical potentials of the 
three-body Hamiltonian ([T]), from which the potential 
that generates the trajectory is subtracted [36| . 

For Coulomb-dominated reactions, the potential that 
generates the classical trajectory is usually the bare P-T 
Coulomb interaction, i.e., 



VpTit, r) = UcT [Rc{t)] + UfT [Rf{t)] 



ZpZtc^ 



(12) 



where Zp and Zp arc the atomic numbers of the projec- 
tile and the target, respectively. 

The TD equation (fTTj) has to be solved for all pos- 
sible trajectories with the boundary condition that the 
projectile is initially in its ground state. 



*'^(t^-oo,6,r) ==0o(r). 



(13) 



This is performed numerically by applying iteratively an 
approximation of the time-evolution operator to the ini- 
tial wave function [35l - l38l |. We use the algorithm de- 
tailed in Refs. [22, [Sa]. At the end of the calculation, 
a breakup probability can be extracted for each trajec- 
tory by projecting the final wave function on the positive- 
energy eigenstates of Hq , 
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Breakup obscrvablcs can be calculated from these prob- 
abilities by proper integration over b [38|. Since these 
observables are obtained by summation over breakup 
probabilities and not over breakup amplitudes, the time- 
dependent technique cannot account for quantum inter- 
ferences between different trajectories. We will see in 
Sec. HV] the effects of such interferences. 



D. Dynamical eikonal approximation 

More recently, the DEA has been developed from the 
comparison between the time dependent model and the 
eikonal approximation [40, |4l[ . It relies on the eikonal 
factorization of the three-body wave function (jB]) [42] 

^^^^{R, r) = e'^"^$(fi, r). (15) 

At sufficiently high energy, the deviation from the initial 
plane wave e'^°^ of the P-T relative motion is expected 
to be small. The dependence on i? of Vf is thus expected 
to be smooth. This enables us to neglect its second-order 
derivative in R with respect to its first-order derivative 

Ajj$(i2, r) < Kod/dZ^{R, r). (16) 

Therefore, introducing the factorization (jlSp into the 
three-body Schrodinger equation ([J), leads to the DEA 
equation [41 1 



^PT OZ 

[{Ho - Eo) + UMRc) + UfTiRf)] $(^, b,r),(17) 

where the dependence of the wave function on the longi- 
tudinal Z and transverse b parts of the projectile-target 
coordinate R has been made explicit (see Fig.[T]). 

The DEA equation ((T7|) is mathematically equiva- 
lent to a time-dependent Schrodinger equation for a 
straight-line trajectory [see Eq. (|TT|)]. It can therefore 
be solved using similar numerical techniques as in the 
time-dependent model |35l - [38{ . As explained in Ref. [4l| , 
we use the algorithm detailed in Ref. [S^. The solution 
is obtained for each transverse component b of the P-T 
coordinate with the boundary condition 
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Breakup amplitudes can then be extracted from the wave 
function 
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Since no semiclassical approximation has been made to 
derive Eq. (|17l) , the coordinates Z and b are quantal vari- 
ables. This enables us to take into account interferences 
between trajectories. This is noticeable by the fact that 
breakup observables are obtained in the DEA by integrat- 
ing amplitudes (|19p over b and not breakup probabilities 
as in the time dependent method of Sec. Ill CI The do- 
main of validity of the approximation (J16p remains to be 
tested by comparison to the other reaction models. 



III. MODEL INPUTS 

In such a comparison, it is important to ensure that the 
inputs of all models are consistent. In this section, we 
summarize the various parameters that have been con- 
sidered in the calculations presented in Sec. IIVI More 
details can be found in the Addendum given as supple- 
mental material. 

All masses are calculated as mass number times the 
nucleon mass toat = 931.5 McV/c^. The effective po- 
tentials simulating the interactions between ^^C, n, and 
Pb are chosen identical in all three models. For the ^^C- 
n potential we take a Woods-Saxon central form factor 
(with depth K;s = 63.023 MeV, radius Rq = 2.651 fm, 
and diffuseness a = 0.6 fm) plus a spin-orbit term (with 
depth Vso = 23.761 MeVfm^, and the same radius and 
diffuseness as the central term). This potential repro- 
duces the two bound states of ^^C: the 1/2+ ground 
state in the lsi/2 partial wave at the experimental en- 
ergy -Eisi/2 = —1.218 MeV and the excited 5/2+ excited 
state as a 0^5/2 state at £^0^5/2 = —0.478 MeV. 

The interactions between the Pb target and the pro- 
jectile constituents are simulated by optical potentials 
chosen from the literature. The Bccchetti and Green- 
lees parametrization [46| is used for the n-Pb poten- 
tial. Since no ^^C-Pb potential is available, we use, at 
68 MeV/nucleon, a potential reproducing the elastic scat- 
tering of ^^O on Pb at 94 MeV/nucleon [ill and, at 
20 MeV/nucleon, a potential fitted to the elastic scatter- 
ing of 1^0 on Pb at 312.6 MeV (potential 13 of Ref. H). 
At both energies, the radius of the potential is scaled 
by (I4I/3 + 208i/3)/(16i/3 + 208^/3) to correct for the 
difference between the sizes of ^^C and ^^O. 

Convergence is an important part of the study and 
therefore was thoroughly tested for all cases. The pa- 
rameter sets quoted in the Addendum provided as sup- 
plemental material ensure an accuracy of at least 4% in 
the energy and angular distributions for all three models. 



IV. RESULTS 

In this section we present the comparison for the 
breakup of ^^C on ^'^^Pb at two different beam energies: 
the first (68 MeV/nucleon) corresponds to an energy typ- 
ical of fragmentation facilities for which data already ex- 
ist [3], and the second (20 MeV/nucleon) serves as an 
example of the energies that will be available in facilities 
such as SPIRAL2 and FRIB. We present the breakup 
cross section as a function of either the '^'^G-n relative 
energy E or the scattering angle of the ^^C-n center of 
mass system. 

The results for ^^C on Pb at 68 MeV/nucleon are pre- 
sented in Figs. 121 and 121 All three models predict nearly 
identical energy distributions (see Fig. [2]) : they differ by 
only 1-3% at the peak. They are also in excellent agree- 
ment with the RIKEN data [J|, validating the reaction 
theory and the assumed single-particle nature of ^^C. The 
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FIG. 2: Distribution for the breakup of ' C on Pb at 
68 MeV/nucleon as a function of the ^*C-n relative energy. 
Comparison of three models: CDCC (solid), TD (dotted), 
and DEA (dashed). Experimental data from Ref. [y|. 
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FIG. 3: Breakup cross section for ^^C on Pb at 
68 MeV/nucleon as a function of the scattering angle of the 
^*C-n center of mass. 
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FIG. 4: Distribution for the breakup of ^'''C on Pb at 
20 MeV/nucleon as a function of the ^*C-n relative energy. 
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FIG. 5: Breakup cross section for '^^C on Pb at 
20 MeV/nucleon as a function of the scattering angle of the 
^*C-n center of mass. In addition to CDCC, TD, and DEA re- 
sults, a TD calculation using straight-line trajectories (dash- 
dotted) is shown. 



aim of this analysis being to compare theories to each 
other, the theoretical cross sections have not been folded 
with the experimental energy resolution. Such a folding 
does not affect much the agreement between theory and 
experiment in the present case. 

The angular distributions are shown in Fig. |31 The two 
quantal models, CDCC and DEA, agree very well with 
each other. In particular they exhibit similar diffraction 
patterns. The TD model does not exhibit any diffractive 
pattern. This diffraction pattern is a quantal effect cor- 
responding to interferences between trajectories, an ef- 
fect excluded in the semiclassical approximation. Never- 
theless, the TD calculation reproduces the general trend 
of the angular distribution at forward angles. This ex- 
plains why, once integrated over the scattering angle, it 
produces a cross section nearly identical to the quantal 
models. Although DEA provides a good approximation 
to CDCC, a slight shift of about 3% in Q is observed be- 
tween both oscillatory patterns. However, at such beam 
energy, this small discrepancy is negligible compared to 
the uncertainties in the optical potentials. 

Next we analyze the breakup of ^^C on Pb at 



20 MeV/nucleon. The energy distribution is displayed 
in Fig. |4] and the angular distribution in Fig. [5] If we 
first focus on the comparison of TD and CDCC models, 
excellent agreement in the energy distribution is found 
(a mere 1% difference at the peak). At this energy too, 
the semiclassical approximation fails at reproducing the 
correct diffraction pattern seen in the CDCC angular 
distribution, but, as seen at higher energy, the general 
trend of the cross section is well approximated by the TD 
model at forward angles. These results show that the TD 
model provides accurate breakup observables integrated 
over the scattering angle even at energies below the range 
of validity mentioned by Alder and Winther [331 ■ Because 
of its semiclassical approximation, the TD model cannot 
account for quantal interferences in the angular distribu- 
tions. Nevertheless, it produces a qualitative estimate of 
the behavior of such distributions at forward angles. 

The DEA energy distribution does not agree with the 
other two models at 20 MeV/nucleon; it is about 10% 
too high at the peak. Due to its quantal nature, DEA 
does exhibit a diffraction pattern in the angular distri- 
bution, but the small discrepancy with CDCC found at 



68 MeV/nuclcon is now significantly increased as the 
DEA angular distribution peaks at more forward angles. 
The shift reaches here 10% in 6. These results suggest 
that the difference observed between DEA and the other 
two models at low energy comes primarily from the lack 
of Coulomb deflection in DEA: Relying on the eikonal ap- 
proximation, the DEA assumes that the incoming plane- 
wave motion of the projectile is not much perturbed by 
its interaction with the target (|16p . The DEA thus forces 
the projectile straight ahead into the high-field zone of 
the target, leading to a larger breakup cross section and 
a more forward angular distribution. On the contrary, 
the usual TD approach, being based on Coulomb tra- 
jectories, naturally includes the Coulomb deflection and 
hence reproduces CDCC calculations fairly well. 

To test this hypothesis, we first repeat the time- 
dependent calculation using straight-line trajectories in- 
stead of hyperbolas (dash-dotted line in Fig. [5]). Of 
course, this TD calculation does not exhibit any diffrac- 
tion pattern. However, it provides a fair approximation of 
the general trend of the DEA angular distribution in the 
same way the usual TD calculation follows the CDCC 
one (compare the dotted and solid lines). This result 
was to be expected as the DEA equation ([T7|) is math- 
ematically equivalent to a time-dependent Schrodinger 
equation with straight-line trajectories (see Sec. IIIDj) . 
It nevertheless confirms the significance of Coulomb de- 
flection in the reaction process. Second, we compare 
DEA to CDCC in a purely nuclear calculation, i.e. set- 
ting Zt = 0. The corresponding angular distributions 
are shown in Fig. [51 At large angles, both calculations 
are nearly identical. At forward angles, however, DEA 
underestimates CDCC and exhibits an oscillatory pat- 
tern shifted to larger angles. This difference with the 
Coulomb-dominated reaction is not very surprising as 
the nuclear interaction, being mostly attractive, tends 
to deflect the projectile within the high-field zone of the 
target. This very stringent test indicates that Coulomb 
deflection is not the only reason for the discrepancies ob- 
served in Figs. ID and [5] and that other effects, such as 
nuclear deflection and/or couplings between various im- 
pact parameters 6, are also signiflcant. 

These results confirm that at low energies, the approx- 
imation ()16p can no longer be performed as it suppresses 
part of the deflection of the projectile by the target, 
and/or sonic coupling effects between different te. A cor- 
rection of the DEA that could account for the Coulomb 
deflection would most likely provide a better description 
of Coulomb-dominated reactions at low-energy. 

To make sure that the qualitative features of our 
analysis do not depend on the particular choice of the 
core-target interaction, we have repeated the Coulomb- 
breakup calculations at 20 MeV/nucleon using the opti- 
cal potentials used at 68 MeV/nucleon. As expected, the 
cross sections are sensitive to the parametrization of these 
potentials. The 68 MeV/nucleon potentials changes the 
energy distribution by 2-5% and reduces the amplitude 
of the oscillations of the diffraction pattern of the angu- 
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FIG. 6: Angular distributions for a hypothetical purely nu- 
clear collision between ^^C on Pb at 20 MeV/nucleon. DEA 
is compared to CDCC. 



lar distribution. It also shifts that pattern by about 3% 
to larger angles. Nevertheless, the qualitative differences 
between the three models remain very similar. 

The difference between CDCC and DEA seems to 
evolve smoothly when reducing the beam energy. At 
40 MeV/nuclcon, the energy distributions differ by a 
mere 5% at the maximum and the shift in the angular 
distributions is also about 5% in 9. 



V. CONCLUSIONS 

In this study, we perform a comparison of non- 
perturbative models of reactions involving loosely-bound 
nuclei. We compare the continuum-discretized coupled 
channel method (CDCC), the time dependent approach 
(TD) and the dynamical eikonal approximation (DEA) 
for the dissociation of a one-neutron halo nucleus on a 
heavy target. Starting from exactly the same three-body 
Hamiltonian, we calculate the energy distribution and 
angular distribution following the breakup of ^^C on Pb 
at 68 MeV/nucleon and 20 MeV/nucleon in all three 
frameworks. 

Our results show that for angle- integrated observables, 
TD works well and can be safely used in the analysis 
of data obtained at both intermediate-energy and low- 
energy facilities, i.e., on an energy range much larger 
than suggested in the original semiclassical approxima- 
tion of Alder and Wither [39|. However, due to its classi- 
cal treatment of trajectories, TD cannot account for the 
diffraction pattern seen in the angular distributions. It 
provides only the general trend of these cross sections at 
forward angles. 

The DEA approach is able to accurately re- 
produce the CDCC angular and energy distribu- 
tions at 68 MeV/nucleon and therefore provides a 
computationally-efficient alternative to CDCC without 
sacrificing accuracy. In contrast, at the lower beam en- 
ergy, both energy and angular distributions in DEA can- 
not reproduce the CDCC results. DEA overestimates the 



energy distribution by 10% and, although the DEA an- 
gular distribution exhibits a diffraction pattern similar 
to that of CDCC, this pattern is shifted to more forward 
angles by about 10%. The primary cause of these dis- 
crepancies is the approximation (|16p made in DEA. For 
Coulomb-dominated reactions, it amounts mostly to the 
absence of Coulomb deflection in that model. Thanks to 
the present analysis we now understand how to remedy 
the problem so that the domain of validity of DEA can 
be partially extended to the lower energies. 

Although valid at all energies, CDCC is a reaction 
model that requires signiflcant computational power. 
Our analysis shows that, depending on the beam energy 
and/or the observable considered, it can be reliably re- 
placed by the TD model or the DEA. Since the DEA and 
TD techniques are computationally less expensive, these 
could allow for improving the description of the projectile 
in reaction models at a reasonable cost. 

The present study corresponds to the first comparison 
of non-perturbative breakup models at intermediate en- 
ergies. It quantitatively shows for which observables and 
energies the models agree and in which conditions their 
predictions should be considered with caution. This pro- 
vides for the first time the range of validity of the three 
models. The projectile description being quite general, 
these results can be extended to other neutral looscl- 
bound systems with confidence. Note that they cannot 



be readily extended to charged systems as the mecha- 
nism of the Coulomb dissociation of proton-halo nuclei 
differs from that of neutron halos: the former involves 
more significant E2 transitions and stron ger higher-order 
dynamical effects than the latter J24| - l26i ISOl Isif . 

Being focused on the comparison between three reac- 
tion models, the present study has been performed within 
the framework of non-relativistic quantum mechanics. 
However, relativistic effects may start to play a role at 
energies around 100 MeV/nucleon [4^. Our conclusions 
should therefore not readily be extended to such energies. 
A detailed analysis of the effect of relativity in breakup 
reactions is planned in the near future. 
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